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Abstract 

In this paper, we introduce a new class of bivariate distributions called the bivariate 
exponentiated extended Weibull distributions. The model introduced here is of Marshall- 
Olkin type. This new class of bivariate distributions contains several bivariate lifetime 
models. Some mathematical properties of the new class of distributions are studied. We 
provide the joint and conditional density functions, the joint cumulative distribution func¬ 
tion and the joint survival function. Special bivariate distributions are investigated in some 
detail. The maximum likelihood estimators are obtained using the EM algorithm. We 
illustrate the usefulness of the new class by means of application to two real data sets. 
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1 Introduction 


The Weibull distribution has assumed a prominen t pos ition as statistical model for data from 


reliability, engineering and biological studies (jMcCooll . 


2 01 2). The Weibull distribution is a 


reasonable choice due to its negatively and positively skewed density shapes. However, this 
distribution is not a good model for describing phenomenon with non-monotone failure rates, 
which can be found on data from applications in reliability studies. Thus, extended forms of 
the Weibull model have been sought in many applied areas. As a solution for this issue, the 
inclusion of additional parameters to a well-defined distribution has been indicated as a good 
methodology for providing more flexible new classes of distributions. 
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The class of extended Weibull (EW) distributions pioneered by Gurvich et al. (j.1997) has 
achieved a prominent position in lifetime models. Its cumulative distribution function (cdf) is 
given by 

G(x; \,£) = l-e~ XH{x *\ x>0, A > 0, (1.1) 

where is a non-negative monotonically increasing function which depends on the pa¬ 

rameter vector The corresponding probability density function (pdf) is given by 


g (x; A, £) = Xh (x; £) e XH ( x ’€\ x > 0, A > 0, 


( 1 . 2 ) 


where h(x;£) is the derivative of H(x;£). We emphasize that several distributions could be 
expressed in the form (11.11). Ta ble [T1 sum marizes several o f these models. Further, we refer the 


reader to 


Nadaraiah and Kota ( 20051 ) and lPham and Lail ( 20071 ) . 


In recent years, many auth ors wor 


ced on this class of distributions such as the beta ex¬ 


tended Weibu 
by 


by. 


Silva et al. 


Cordeiro et al 


(2012), the extended Weibull power series distributions 


1 family by ( 

(2013), the complementary extended Weibull power series class of distributions 


Cordeiro and Si 


Santos-Neto et al. 


butions by 


val feOTiJ), the Marshall-Olkin extended Weibull family of distributions by 


(2014) and the expon entiated extended Weibull-power series class of distri- 


Tah maseb i and Jafari (12 0151 ). 


The aim of this paper is to introduce a new bivariate exponentiated extended Weibull 
(BEEW) family of distributions, whose marginals are exponentiated extended Weibull (EEW) 
distributions. It is obtained using a meth od s imilar to that used to obtain Marshall-Olkin 


bivariate exponential model (Marshall and Olkm 


1967). The proposed BEEW class of dis¬ 


tributions is constructed from three independent EEW distributions using a maximization 
process. Creating a bivariate distribution with given marginals using this technique is noth¬ 
ing new. The joint cdf can be expressed as a mixture of an absolutely continuous cdf and 
a singular cdf. The joint pdf of the BEEW distributions can take different shapes and the 
cdf can be expressed in a compact form. The joint cdf, the joint pdf and the joint sur¬ 
vival function (sf) are in closed forms, which make it convenient to use in practice. The new 


(Kundu and Gupta, 

2009 

), bivariate generalized linear failure rate ( 

Sarhan et al., 

2011) 


variate generalized Gompertz flE l-Shcrpi env et all l2013l ). bivariate exponentiated generalized 


Weibull-Gom pertz (El-Bassioun et al. 
sion (El-Goharv and El-Morshedv, 2015 


20151 ). bivariate exponentiated modified Weibull exten- 


distributions. This class defines at least 46 (2 x 23) 


bivariate sub-models as special cases. 



























































Table 1: Special cases of EW distributions and corresponding H(x',£) function 


Distribution 

Support 


A 

£ 

Reference 

Exponential 

x > 0 

X 

A 

0 

Johnson et al. (1995) 

Pareto 

x = k 

log {x/k) 

A 

k 

Johnson et al. (1995) 

Gompertz 

x > 0 

c -1 [exp ( cx ) — 1 ] 

A 

c 

GoniDertz 11825) 

Weibull 

x > 0 

X 7 

A 

7 

Frechet (1927) 

Frechet 

x > 0 

X -7 

A 

7 

Frechet 11927) 

Lomax 

x > 0 

log(l T x) 

A 

0 

Lomax [1954) 

Weibull Kies 

0 < n < x < a 

(x - n) bl /(a - x) b2 

A 

(fi,a,bi,b 2 ) 

Kies 11958) 

Log-logistic 

x > 0 

log(l + x c ) 

A 

c 

Fisk 11961) 

Linear failure rate 

x > 0 

ax + bx 2 /2 

1 

(a,b) 

Barlow 11968) 

Log-Weibull 

—OO < X < oo 

exp[(x - n)/a\ 

1 

(M, A) 

White 11969) 

Exponential power 

x > 0 

exp((c*)“ — 1 ) 

1 

(a,c) 

Smith and Bain 11975) 

Burr XII 

x > 0 

log(l +x c ) 

A 

c 

Rodriguez (1977) 

Rayleigh 

x > 0 

x 2 

A 

0 

Ravleieh (1880) 

Phani 

0 < H < x < a 

/ *r 

i 

i 

cr 

A 

(fi, a, b) 

Phani 11987) 

Additive Weibull 

x > 0 

(x/far + (x/fcr 2 

1 

(ai, 02, /3i, /I 2 ) 

Xie and Lai 11995) 

Chen 

x > 0 

expIaJ — 1 ) 

A 

b 

Chen 12000) 

Pham 

x > 0 

{a x f - 1 

1 

(a, P) 

£ham ( 2002 ) 

Weibull extension 

x > 0 

c [expfca :) 6 — l] 

A 

( 7 , b, c) 

Xie et al. (2002) 

Modified Weibull 

x > 0 

a; 7 exp(ca;) 

A 

(7 1 c ) 

Lai et al. (2003) 

Traditional Wibull 

x > 0 

x d exp(cx a — 1 ) 

A 

(a, b, c) 

Nadaraiah and Kotz 12005) 

Generalized Weibull power 

x > 0 

[1 + (x/a) b T - 1 

1 

{a,b,c ) 

Nikulin and Haehiahi (2006) 

Flexible Weibull extension 

x > 0 

exp (ai* — Pi/x) 

1 

{ 011 , Pi) 

Bebbington et al. (2007) 

Almalki Additive Weibull 

x > 0 

ax e + bx^e cx 

1 

( a > b, c, 9, 7 ) 

Almalki and Yuan [2013) 


The usual maximum likelihood estimators can be obtained by solving non-linear equations 
in at least five unknowns directly, which is not a trivial issue. To avoid difficult computation 
we treat this problem as a missing value problem and use the EM algorithm, which can be 
implemented more conveniently than the direct maximization process. Another advantage of 
the EM algorithm is that it can be used to obtain the observed Fisher information matrix, which 
is helpful for constructing the asymptotic confidence intervals for the parameters. Alternatively, 
it is possible to obtain approximate maximum likelihood estimators by estimating the marginals 


first anc 


by (Joe, 


then estimating the dependence parameter through a copula function, as suggested 


1997, Chapter 10), which has the same rate of convergence as the maximum likelihood 


estimators. This is computationally less involved compared to the MLE calculations. This 
approach is not pursued here. Although in this paper we mainly discuss the BEEW, many of 
our results can be easily extended to the multivariate case. 

The main reasons for introducing this new class of bivariate distributions are: (i) This 
class of distributions is an important model that can be used in a variety of problems in 
modeling bivariate lifetime data, (ii) It provides a reasonable parametric fit to skewed bivariate 



















































data that cannot be properly fitted by other distributions, (iii) The joint cdf and joint pdf 
should preferably have a closed form representation; at least numerical evaluation should be 
possible, (v) This class contains several special bivariate models because of the general class of 
Weibull distributions and the fact that the current generalization provides means of its bivariate 
continuous extension to still more complex situations; therefore it can be applied in modeling 
bivariate lifetime data. 

The rest of the paper is organized as follows. We define the EEW and the BEEW class of 
distributions in Section [2l Different properties of this family are discussed in this section. The 
special cases of the BEEW model are considered in Section [3l The EM algorithm to compute 
the MLEs of the unknown parameters is provided in Section U The analysis of two real data 
sets are provided in Section [5j Finally, we conclude the paper in Section [6l 


The BEEW model 


In this s ecti on, we introduce the BEEW distributions using a method similar to that which was 


used by 


Marshall and Olkin (1967) to define the Marshall Olkin bivariate exponential (MOBE) 


distribution. 

First, consider the univariate EEW class of distributions with cdf given by 
Feew (x‘, a, A, £) = ^1 — e~ XH ( x '^ , x > 0, a > 0, A > 0. 
The corresponding pdf is 

/eew {x\ a, A, £) = a\ h (x; £) e - XH ( x & - e -Atf(x;£)^ 


a—1 


( 2 . 1 ) 


( 2 . 2 ) 


From now on a EEW class of distributions with the shape parameter a, the scale parameter 
A and parameter vector £ will be denoted by EEW(a, A, £). Note that many well-known models 


generalizec 

extension 


Rayleigh 


Xie et al 


Surles and Padgett! . 


2 OO 2 L generalized mo dified Weibull dCarrasco et al 


2001 


1993 

), generalized exponential ( 

=»U —- -± 

Gupta and Kundu. 

1999 

), Weibul 

extension ( 

Chen, 

2000) 


Kundu and Raqab, 


2001 


linear failure rate (Sarhan and Kundu 


2009h . generalized Gompert z (lEl-Goharv et al 


and exponentiated modified Weibull extension Sarhan and Apaloo, 


, modified Weibull 


2008 ) generalize d 


2013 ) 


2013) distributions. 


When a is a positive integer, the EEW model can be interpreted as the lifetime distribution 


of a parallel system consisting of a independent and identical units whose lifetime follows the 
EEW distributions. 

































































From now on unless otherwise mentioned, it is assumed that a\ > 0; 02 > 0; 03 > 0 
and A > 0. Suppose U\ ~ EEW(cki,A,£), U 2 ~ EEW(a 2 ,A,£) and C /3 ~ EEW(a 3 ,A,£) and 
they are mutually independent. Here means follows or has the distribution. Now define 
X\ = maxjl/i, C/ 3 } and X 2 = max{[/ 2 , C/ 3 }. Then, we say that the bivariate vector (A}, X2) has 
a bivariate exponentiated extended Weibull distribution with the shape parameters a\, a 2 and 
03 , the scale parameter A and parameter vector £. We will denote it by BEEW(ai, ot 2 , a 3 , A, £). 


Before providing t 
According to 


re joint cdf or pdf, we firs t mention how it may occur in practice. 


Kundu and Gupta (2009), suppose a system has two components and it 


is assumed that each component has been maintained independently and also there is an 
overall maintenance. Due to component maintenance, suppose the lifetime of the individual 
component is increased by Ui amount and because of the overall maintenance, the lifetime 
of each component is increased by U 3 amount. Therefore, the increased lifetimes of the two 
component are X\ = max{C/i,C/ 3 } and X2 = max{t/ 2 , C/ 3 }, respectively. 

We now study the joint cdf of the bivariate random vector (X\ , X2) in the following theo¬ 
rem. 


Theorem 2 . 1 . //(A}, A}) ~BEEW(«i, 02, « 3 , A, £), then the joint cdf of (X\, X2) for x\ > 0 , 
x'2 > 0 , is 

F B EEW (z„z 2 ) = (1 - (j _ e-«<*z «))" (1 - ( 2 . 3 ) 

where z = minjxi, X 2 }. 

Proof. Since the joint cdf of the random variables X\ and X2 is defined as 


FbeewOd, X2) = P (X 1 < xi, X 2 < X2) 

= P (max (Dr, U 3 } < x±, max {C/2, C/3} < X2) 

= P{Ui < xi, U 2 < X2, U 3 < min(xi, X2)). 

As the random variables Ui, (i = 1,2,3) are mutually independent, we directly obtain 


Tbeew(ai, X2\ ai, « 2 , a-3, A, £) = FeewC^i; «i, A, ^)Feew (^25 « 2 , A, ^)Feew(^; 03 , A, ^). (2.4) 
Substituting from 12.11 into 12.41 we obtain 12.31 which completes the proof of the theorem. □ 


Corollary 2 . 1 . The joint cdf the BEEW(ai, 0 : 2 , < 23 , A, £) can also written as 


Freew (#1, X2 ) 


(l_ e -Atf(*i;€))“ 1+ “3(l_ e -A if X1 < X2 

(1 - e ~ XH ^) ai (1 - e-^(^i«) aa+ “3 if Xl > X2 


(2.5) 





-Frew (xi; ai + 03, A, £) Eeew (x 2 ‘, « 2 , A, £) 
TeEW (x\\Oil, A, £) -FeEW (x 2 ‘, «2 + « 3 , A, £) 
-Peew (x; «i + «2 + a 3 , A, £) 


if 

if 

if 


X\ < X 2 
X 2 < Xi 


Xl = X 2 = X. 


The following theorem gives the joint pdf of the random variables X\ and X 2 which is the 


joint pdf of BEEW(«i, « 2 , a 3 , A, £). 


Theorem 2.2. If (X\, X 2 ) ~ BEEW(ai, 02 , a 3 , A, £) then the joint pdf of (X\, X 2 ) for x\ > 0, 


X 2 > 0, is 


where 


/beew (xi, X2) 


/i(xi,x 2 ) if 
< f2(x!,x 2 ) if 
. /o(x) if 


0 < Xl < X 2 
0 < X 2 < Xl 

0 < Xi = X2 = X, 


( 2 . 6 ) 


fi (xi,x 2 ) 


f 2 (xi,X 2 ) 


/o (a:) 


/eew (xi; ai + a 3 , A, £) /eew (x 2 l ol2, A, £) 

(ai + a 3 ) a 2 A 2 h (xi; £) h (x 2 ; 0 (l - e~ XH ^ xv ’^ 

x _ e -XH(x2-4)y 2 ~ 1 e -XH(x 1 -4)-XH(x 2 -4) 

/eew {xi, ai, A, £) /eew (x 2 ; a 2 + a 3 , A, £) 

{ 0 L 2 + a 3 ) aiA 2 h (xi; £) h (x 2 ; £) (l - e _AH(3;i; ^^ 
x _ e -Afr(x 2 ;€)^ aa+ “ 3-1 e -AH(*i;0-AH( a! 2;€) 

0^3 

-;—-—;-/eew (a:; ai+a2 + 03, A, £) 

ai + 02 + a 3 


(2.7) 


( 2 . 8 ) 


(2.9) 


Proof. First assume that xi < x 2 - Then, the expression for f\ (xi,x 2 ) can be obtained simply 
by differentiating the joint cdf Ebeew (xi, X2) given in (12.511 with respect to xi and X2. Similarly, 
we find the expression of f 2 {xi,x 2 ) when x 2 < x\. But /o (x) cannot be derived in the same 
way. Using the facts that 


T OO pX2 TOO (*X\ TOO 

/ / /i(xi,x 2 ) dx\dx 2 + / / /2(xi,x 2 ) dx 2 dx\ + / / 0 (x) dx 

J 0 do do do do 


= 1 , 


/o ao 

(*00 rX 2 


[ X2 Z" 00 / «l+«2+«3 —1 , 

J fi(xi,x 2 ) dx\dx 2 = a-2 J A/i(x;£)^l — e e dx, 


and 


TOO Tail TOO / . 

J J / 2 (xi,x 2 ) dx 2 dxi = ai j A h (x; £) (l - 


ai+Q2+«3 — 1 


- A ZZ(x ;0 


/o ao 
Note that 


/■oo roo 

J fo (x) dx = a 3 J A h (x; £) (l - e ~ XH( - x ’^J 


«1+Cf2+Q!3 —1 


e ~ XH ( x ;£) dx = -—-. 

ai + 0 L 2 + o 3 






Thus, the result follows. 


□ 


The following theorem gives the marginal pdf’s of X\ and X 2 . 

Theorem 2 . 3 . The marginal distributions of X\ and X-2 are EEW(ai + 03, A, £) and EEW (02 
+03, A, ^), respectively. 

Proof. The marginal cdf for X,; is 


Fxi ( Xi ) = P (Xi <Xi) = P (max {Ui, U 3 } < Xi) = P(Ui < X{, U 3 < Xi). 
Since the random variables t/*, (i = 1,2) are mutually independent, we obtain 

Fxi (xi) = P(Ui< Xi) P(U 3 < Xi) 

= -^eew {xi',ai,\,£) Teew (xj;a 3 ,A,£) 

= -^EEW (%i\Oii + 03, A, £) . 

From 12.101 we can derive the pdf of Xi by differentiation. 


( 2 . 10 ) 


□ 


The BEEW model has both an absolute continuous part and a singular part, similar to 
Marshall and Olkin’s bivariate exponential model. The joint cdf of X\ and X2 has a singular 
part along the line x± = x 2 , with weight , u+ q!) + q ., 1 and has an absolutely continuous part on 
0 < X! ^ x 2 < 00 with weight a °^+a 3 • 

Interestingly, the BEEW model can be obtained by using the MarshalLOlkin (MO) copula 
with the marginals as the EEW distributions. To every bivariate cdf Fx i ,x 2 with continuous 
marginals Fx x and Fx 2 there corresponds a unique bivariate cdf with uniform margins C : 
[0, l] 2 —> [0,1] called a copula, such that Fx lt x 2 (a?i, a?2) = C{Fx 1 (xi) ,Fx 2 (^2)} holds for all 
(xi, X2) G M 2 ( Nels on. 1999 ). The MO copula is 


C e i,6» 2 {u ll u 2 ) = u\ 01 u\ 02 min ju® 1 , u 0 2 2 j , 


for 0 < 61 < 1 and 0 < 62 < 1 - Using Ui = Fx t {xf) where Xi is EEW(cq + ct 3 ,A,£) and 
6 i = a ^ a3 , i = 1, 2,3, gives the same joint cdf Fx 1 ,x 2 as m - 

The following result will provide explicitly the absolute continuous part and the singular 
part of the BEEW cdf. 

Theorem 2 . 4 . If {X\. X 2 ) ~ BEEW(ai, 02,03, A, £), then 

Fx i,x 2 (xi,x 2 ) = — Ql + °! 2 — F a (xi, x 2 ) 4- , a3 , - F s (xi,x 2 ), 


O 1 + 0.2 + 0 3 


0\ + o 2 + o 3 












where for x = min {xi, X 2 }, 


F s (xi,x 2 ) 




01 + 0 : 2+03 


and 


F a (xi,x 2 ) = 


Oil + <+2 + O3 


Oi + <+2 
«3 


_ g-Afffxi;^® 1 ^ _ e -A H(x 2 4)y 2 ^ _ e -\H{x£)^ 

(l - e -A"(*;€)) 


01 + 02+03 


Q+ + <+2 

/iere F s (.,.) and F a (.,.) are t/ie singular and the absolute continuous parts, respectively. 


Proof. To find F a (xi,x 2 ) from F Xl ,x 2 (xi,x 2 ) = aF a (xi, x 2 ) + (1 - a) F s (xi, x 2 ) , 0 < a < 1, 
we compute 

d 2 F Xl ,x 2 {xi,x 2 ) I fi(x l 7 x 2 ) if xi < x 2 

ox 1 dx 2 [ f2(xi,x 2 ) if Xi > x 2 , 

from which a may be obtained as 

f OO [-X2 


a = 


f°° f X2 f Xl ai + a 2 

/ / fi(xi,x 2 ) dx\dx 2 + / / f 2 (xi,x 2 ) dx 2 dxi = 

Jo Jo . 


and 


/o Jo 


n X2 

/a (s, f) ds dt. 


Ot 1 + (+2 + (+3 


Once a and F a (.,.) are determined, F s (.,.) can be obtained by subtraction. 


□ 


Corollary 2 . 2 . The joint pdf of X\ and X 2 can be written as follows for x = min {xi, x 2 }- 
fx i,X 2 (xi,X 2 ) = ——^——^—fa(x 1,X 2 ) + —— —fs ( X ) 


01 + a 2 + 03 


<+1 + <+2 + 03 


where 


j. ^ x ^ _ a i + a 2 + a 3 ^ j /eew (•'Ci; «i + 03, A, £) /eew (+ 2 ; <+ 2 , A, £) if xi < x 2 , 
«i + (+2 1 /eew (xi', 011 , A,£) /eew (+2; oi 2 + 03 , A,£) if xi > X 2 , 


and 


/s (x) = /eew (+ 01+02+ 03, A, £) . 


Clearly, here f a (xi,x 2 ) and f s (x) are the absolute continuous part and singular part, respec¬ 
tively. 


Having obtained the marginal pdf of X\ and X 2 , we can now derive the pdf’s as presented 


in the following theorem. 



Theorem 2 . 5 . The conditional pdf of X % given Xj = Xj, denoted by fx^x, {xi\xj), i ^ j = 1,2, 
is given by 


fXi\Xj (Xi\Xj) = < 


fxi\Xj (■ x i\ x 3 ) if 0 < 

/£W, if 0 < x 3 < X i 

fx]\ Xj ( x i\ x i ) if x i = x i > °> 


( 2 . 11 ) 


where 




OL 3 


(i) (a* + a 3 )<x,A h(xi;$,) (l-e 

J Xi\ Xj W X i> ~ ( a?2 + a3 ) (1 - 

fx-\ Xj i x i\ x j ) = h ( x f,£) (i - ‘ e~ XH ( Xi ’® 

3) fa.U..)- a 3 / x -Agfe;£)\ Qi 

- aj - + a 3 V J ' 

Proof. The proof follows readily upon substituting the joint pdf of (X\, X 2 ) given in Theorem 
12.21 and the marginal pdf of Xj, given in Theorem 12.31 using the following relation 


t ( | fx i)Xj {Xi,Xj) . 

fx,lx, fefcf) = {xi) . * = 1.2. 


( 2 . 12 ) 

□ 


Proposition 1. Since the joint sf and the joint cdf have the following relation 


S Xl ,x 2 (xi,x 2 ) = 1 - F x 1 (xi) - Fx 2 (x 2 ) + F Xl ,x 2 (xi,x 2 ) , 


(2.13) 


therefore, the joint sf of X\ and X 2 also can be expressed in a compact form. 


Proposition 2 . 


Basu tl 97 n) defined the bivariate failure rate function hx 1 ,x 2 (x\,x 2 ) for the 


random vector (X 3 ,X 2 ) as the following relation 


hx!,x 2 (xi,x 2 ) = 


fx!,x 2 (xi,x 2 ) 


(2.14) 


S Xl ,x 2 (xi,x 2 ) 

We can obtained the bivariate failure rate function hx 1 ,x 2 (xi,x 2 ) for the random vector (X\, 
X 2 ) by substituting from (12.61) and (12.131) in (12 . 14[) . 


Lemma 2 . 1 . The cdf ofY = niax{Xi, X 2 \ is given as 


F y (y) = (l - e ~ XH ^) 


ai+a 2 +a3 


(2.15) 


Proof. Since 


F Y (y) = P(Y<y) = P(max{X 1 ,X 2 }<y) 


P {Xi <y,X 2 <y) = P (max{f/i, U 3 } < y, max{[/ 2 . U 3 } < y) 














p {U\ <y,U 2 < y,U 3 < y ), 


where the random variables Ui (i = 1,2,3) are mutually independent, we directly obtain the 
result. □ 

Lemma 2 . 2 . The cdf of T = min{Xi, Al 2 } is given as 

F T ( t ) = (l - e -^))“ 1+ “ 3 + (l - e -^))“ 2+ " 3 _ (i _ e ~ mwy ^ +a \ 

Proof. It is easy to prove that by using Equations (12.13|) and (12.151) . □ 


3 Special cases 

In this Section, we consider some special cases of the BEEW distributions. 


3.1 Bivariate generalized exponential distribution 

If H (x;£) = x, then the joint cdf (12.5)) becomes 

r (i-e-r^u-e-^r * *i <* 2 

1 15 2J 1 (l- e -^) ai (l-e - A ^)“ 2+ “ 3 if x 1>X2 , 


which is the joint 


cdf of bivariate generalized exponential (BGE) distribution introduced by 


Kundu and Gupta (2009). By Theorem 12.101 the marginal distributions of X\ and X 2 are 


GE(aq + 03 , A) and GE (a 2 + 03 , A), respectively. 


3.2 Bivariate generalized linear failure rate distribution 

If H ( x ; £) = fix + 3) x 2 and A = 1, then the joint cdf (12.51) becomes 

(l - e -^~ 2 x iy i+a3 (l - if Xl <x 2 

F{j ' ] " r2) { ^ 1 - e -^* 1 -? a! ?) a i(l-e-^* a -? x a) a2+ “3 if xi>x 2 , 


which is the joint cdf o: 


duced by 


Sarhan et al. 


biv ariate generalized linear failure rate (BGLFR) distribution intro- 


(|201ll ). By Theorem 12. 101 the marginal distributions of X\ and X 2 are 


GLFR(ai + 03 , (3, 7 ) and GLFR (a 2 + 03 , /3, 7 ), respectively. 


3.3 Bivariate exponentiated Weibull distribution 

If H (x;£) =x@, then the joint cdf (12. 5 p becomes 


F (xi, x 2 ) = 


l—e 


1—e 


-A*?j“ 1+a3 ^ 1 _ e -^ aa 

-A,?rvi_ e -Ax?r 2+as 


) 


if x'i < x 2 
if xq>x 2 . 
















We call this, bivariate exponentiated Weibull (BEW) distribution. By Theorem 12.101 the 
marginal distributions of X\ and X 2 are EW(au+a 3 , A, f3) and EW («2 + « 3 , A, f5) , respectively. 


3.4 Bivariate generalized Gompertz distribution 

If if(x ;0 = r 1 ( e@ x — 1), then the joint cdf (12.51) becomes 

f (l _ e -A^ 1 (e^-l)) ai+ “ 3 (l_ e -A/3- 1 (e^-l ) r 2 if Xl < X2 

F(X " " 2 *”\ if Xl> ^ 

which is th e jo int cdf of bivariate generalized Gompertz (BGG) distribution introduced by 


El-Sherpienv et al. 


(201^). By Theorem 12.101 the marginal distributions of X\ and X 2 are 


GG(ai + 03 , A, (3 ) and GG (02 + & 3 , A, /3 ), respectively. 


3.5 Bivariate exponentiated generalized Weibull-Gompertz distribution 


If H (x;£) = x^(e yx5 — 1), then the joint cdf (12.51) becomes 
F (x'i, x 2 ) = < 


S \ ai+a3 / 0 6 \ a2 

1 _ g-Ax^e 7 i-i) j ( X _ e -Ax 2 (e 7 2-1) | j£ X\ < X2 


S \ a l / a s \ “ 2+«3 

1 _ g-Axj’le 7 i-i) j ( £ _ e -Ai 2 (e 7 2-1) | )f Xl > X2j 


which is the joint cdf of bivariate exponentiated generalized Weibull-Gompertz (BEGWG) dis¬ 


tribution introduced by 


El-Bassioun et al. 


(2015). By Theorem l2.1Q[ the marginal distributions 


of Xi and X 2 are EGWG(ou + 03 , A, /3, 7 , 5) and EGWG (02 + a 3 , A, j3, 7 , 5) , respectively. 


3.6 Bivariate exponentiated modified Weibull extension distribution 

If H (x;£) = /3(e( x /PF — 1), then the joint cdf (12.51) becomes 

if Xl < X2 
if XI>X2i 


which is the joint cdf of bivariate exponentiated moc 


tribution introduced by 


El-Goharv and El-Morshedv 


ihed Weibull extension (BEMWE) dis- 


2015). By Theorem 12.101 the marginal 


distributions of X\ and X 2 are EMWE(ai + 03 , A, f3, 7 ) and EMWE (a .2 + 03 , A, /3, 7 ), respec¬ 
tively. 


4 Maximum likelihood estimation 

In this section, we first study the maximum likelihood estimations (MLE’s) of the parameters. 
Then, we propose an Expectation-Maximization (EM) algorithm to estimate the parameters. 























Let (xn,xi 2 ), • • •, {xi n ,x 2n ) be an observed sample with size n from BEEW distribution 
with parameters 0 = (ctu, a 2 , a 3 , A, £) . Also, consider 


h = {i : xu = x 2i = Xi} , h = {i: x u < x 2i }, h = {i ■ x u > x 2 i }, i = 1 , 

and 

n 0 = |/o|, ni = |/i|, n 2 = \I 2 \, n = n 0 + ni + n 2 . 
Therefore, the log-likelihood function can be written as 

us) = Xi log (/l i x li, X 2 i)) + 22 log (f 2 (xii, X 2 i)) + ^ log (/o (xj)) 
iel 1 iS/2 iG/o 

= (2m + 2n 2 + no) log (A) + mlog (a 2 ) + n 2 log (aq) + n 0 log (a 3 ) 

+nilog(ai+ a 3 )+ n 2 log(a 2 + a 3 ) + 22 log {h (xh; £)) 

ie/iu /2 

+ 22 !°g + '22^0g( h (x i -^)) 

iEliUl 2 iG/o 

+ (ou + a 3 — 1) I Y, log (l - e~ XH ^ Xu '^ + 22 lo S (l- _ e - XH ( X2i &^ 
\ieh ieh 

+ (a 2 - 1) 22 log (l - e~ AH ( X2i; &) + (aq - 1) 22 log (1 _ e ~ XH ( Xli & 


,n, 


ieh 


ieh 


+ (aq + a 2 + a 3 - 1) 22 lo g [} ~ e 
ieh 

+A E Xi + ^ xi i+ 22 X2i 

VzG/o iehul 2 2 G-/ 1 U /2 


(4.1) 


where / 1 , f 2 and /o are given in (12.71) . (12.81) and (12.91) . respectively. We can obtain the MLE’s 
of the parameters by maximizing £ (0) in (14.11) with respect to the unknown parameters. This 
is clearly a (k + 4)-dimensional problem. However, no explicit expressions are available for the 
MLE’s. We need to solve (k + 4) non-linear equations simultaneously, which may not be very 
simple. Therefore, we present an expectation-maximization (EM) algorithm to find the MLE’s 
of parameters. It may be noted that if instead of (Aq,X 2 ), we observe U\, U 2 , and I/ 3 , the 
MLE’s of the parameters can be obtained by solving a two-dimensional optimization process, 
which is clearly much more convenient than solving a (fc + 4)-dimensional optimization process. 
For this reason, we treat this problem as a missing value problem. 

Assumed that for the bivariate random vector (Aq,X 2 ), there is an associated random 
vectors 


Ai = 


0 U\ > U 3 
1 t/i < t / 3 


and A 2 = 


0 U 2 > U 3 
1 U 2 <U 3 . 


Note that if X\ = X 2 , then Ai = A 2 = 0. But if X\ < X 2 or X\ > X 2 , then (Ai,A 2 ) 
is missing. If (Xi,X 2 ) € I\ then the possible values of (Ai,A 2 ) are (1,0) or (1,1), and If 
(Xi, X 2 ) € I 2 then the possible values of (Ai, A 2 ) are ( 0 , 1 ) or ( 1 , 1 ) with non-zero probabilities. 

Now, we are in a position to provide the EM algorithm. In the E-step of the EM-algorithm, 
we treat it as complete observation when they belong to Iq. If the observation belong to I\, 
we form the ‘pseudo’ log-likelihood function by fractioning (x\,X 2 ) to two partially complete 
“pseudo” observations of the form (x\,X 2 ,ui ( 0 )) and (xi, a; 2 , it 2 ( 0 )), where u\ ( 0 ) and u 2 ( 0 ) 
are the conditional probabilities that (Ai,A 2 ) takes values (1,0) and (1,1), respectively. It is 
clear that 


Ul ( 0 ) 


o 11 

OL\ + Q!3 


u 2 (0) 


«3 

«1 + Q!3 


Similarly, If the observation belong to / 2 , we form the ‘pseudo’ log-likelihood function 
of the from (yi,y 2 ,ui (0)) and (xi,x 2 ,u 2 (0)), where v\ (0) and u 2 (0) are the conditional 
probabilities that (Ai,A 2 ) takes values (0,1) and (1,1), respectively. Therefore, 


vi ( 0 ) 


0.2 

0-2 + «3 ' 


u 2 (0) 


«3 

Ot 2 + «3 


For brevity, we write u\ (0), tt 2 (0), v\ (0), u 2 (0) as ui, u 2 , v\, u 2 , respectively. 


E-step: Consider hi = E(N\yn, y 2 j, 0). The log-likelihood function without the additive 
constant can be written as follows: 


^pseudo (©) 


(no + 2 ni + 2n 2 ) log (A) + (u\n\ + n 2 ) log (cci) + (n\ + v\ n 2 ) log (cc 2 ) 

+ (n 0 + u 2 n 1 + v 2 n 2 ) log (0:3) + X log (h (xf, £)) + Y log (h (x u \ £)) 

i£lo iGl\Ul 2 

+ X l °s(h(x2u^)) + (ai + a2 + a 3 -l)Y lo s{ 1 -e~ XH{xi ’^ 

iehLIh ielo 

+ (ai + <33 - 1) X l°g (i — 
ieh 

+ (a 2 + a 3 - 1) X log (l - e~ XH ( X2i ’^ 
ieh 

+ (a 2 - 1) X log ^1 - + (ai - 1) X log (l - e - XH ^'^ 

i£l 1 26/2 


-A 


X]- H ’( x *^) + X! X ^( x 2i;C) 


1 ie/o 


ie/iU/2 


iehuh 


M-step: At this step, ^ pse udo (0) is maximized with respect to ai,a 2 ,a 3 , A and For fixed 



A and £, the maximization occurs at 


ai(A,0 


a 2 (A) 


a 3 (A) 


uini + n 2 

Y.,a, w ^) + Y.,a t ui,W(xK)' 

Til + VlU 2 

n 0 + u 2 ni + u 2 n 2 

Eie/o ^(*0 + Eie/i + E j6 / 2 ^(X2i) ’ 


(4.2) 

(4.3) 

(4.4) 


where W(x) = log (l — e XH ( x &y For fixed 01 , 0 : 2 , 0:3 and £, ^ pse udo (©) is maximized with 
respect to A as a solution of the following equation: 


where 


n 0 + 2 ni + 2 n 2 
g(A) 


A, 


(4.5) 


_ , , , ( l4) 

g(A) - -(oi +0 2 + 0 3 - l)^ X _ Q—\H(xi\£,) 

ieio 

- (O! + «3 - 1) 2. ! _ - ( " 2 + «3 - 1) 2^ ! _ e -\H(x 2 ^) 

id 1 iE /2 

, ^^H(x 2 i ;Z)e- XH (*^ H ( Xl ..£) e -mxu-,& 

- (“2 - X ) Z, l_ e -Ag (3;2< ;0 - (“1 - !) Z, 

ie/i ie+ 

+ ^ff(xj;^)+ ^ #(a+b£) + ^ H{x 2 i \i). 

id 0 iE/iU /2 ie/iU /2 


Finally, for fixed oi, o 2 ,03 and A, ^ pseu do (©) is maximized with respect to £ as a solution 
of the following equation: 

^ £ pseudo (©) = 0. (4.6) 

The following steps can be used to compute the MLE’s of the parameters via the EM 
algorithm: 

Step 1: Take some initial value of ©, say 0^ = ^oj°\ oj>°\ o^, A^, ^ 0) j . 

Step 2: Compute u\, u 2 , v \, and v 2 . 

Step 3: Find A by solving the equation (14.51) . say A^k 
Step 4: Find £ by solving the equation (14.61) . say ^ \ 

Step 5: Compute 0 ^ = dj(A^\£^), * = 1,2,3 from (I4.2I) - (I4.4[) . 

Step 6: Replace 0^°^ by ^), go back to step 1 and continue the 

process until convergence take place. 











5 Two real examples 


We consider BEEW distributions for fitting these two data sets. But, this family of distributions 
is a large class of distributions. Here, we consider six sub-models of BEEW distributions: BGE, 
BGLFR, BEW, BGG, BEGWG, and BEMWE. Some of them are suggested in literature. 

Using the proposed EM algorithm, these models are fitted to the bivariate data set, and 
the MLE’s and their corresponding log-likelihood values are calculated. The standard errors 
(s.e.) based on the observed information matrix are obtained. 

For each fitted model, the Akaike Information Criterion (AIC), the corrected Akaike in¬ 
formation criterion (AICC) and the Bayesian information criterion (BIG) are calculated. We 
also obtain the Kolmogorov-Smirnov (K-S) distances between the fitted distribution, the em¬ 
pirical distribution function, and the corresponding p-values (in parenthesis) for X\, X 2 and 
max(Xi, X 2 ). Finally, we make use the likelihood ratio test (LRT) and the corresponding 
p-values for testing the BGE against other models. 


Meintanis (2007) and is obtained from the group 


Example 5.1. The data set is given from 
stage of the UEFA Cha mpio n’s League for the y ears 2004-05 and 2005-2006. In addition, 


Kundu and Gupta 12000) and 


Sarhan et al. 


201A ) analyzed this data. The data represent the 


football (soccer) data where at least one goal scored by the home team and at least one goal 
scored directly from a kick goal (like penalty kick, foul kick or any other direct kick) by any 
team have been considered. Here X\ represents the time in minutes of the first kick goal scored 
by any team and X 2 represents the first goal of any type scored by the home team. 

The results are given in Table [ B It can be concluded that all six models are appropriate for 
this data set. But, the BGW and BGG distributions are better than other distributions. 


Example 5.2. The dat a set was first published in ‘Washington Post’ and is available in 


Gsorad and Welsh 1989) . It is represent the American Football League for the matches on 


three consecutive weekends in 1986. Here, X\ represents the ‘game time’ to the first points 
scored by kicking the ball between goal posts , and represents the ‘game time’ to the first points 


scored , by m oving the ball into the end zon e. \Kundu and Gupta 201(1) 


Jamalizadeh and Kundu 


(2013). and 


Balakrishna and Shin (2014) analyzed this data. We divided all the data by 100. 


The results are given in Tabled It can be concluded that all six models are appropriate for this 
data set. But, the BGE distribution is better than other distributions. 























Table 2: The MLE’s, log-likelihood, AIC, AICC, BIC, K-S, and LRT statistics for six sub¬ 
models of BEEW distribution of first data set. 























































































Table 3: The MLE’s, log-likelihood, AIC, AICC, BIC, K-S, and LRT statistics for six sub¬ 
models of BEEW distribution of second data set. 



Distribution 

Statistic 

BGE 

BGLFR 

BEW 

BGG 

BEWG 

BEMWE 

bi 

0.0921 

0.0921 

0.1367 

0.0921 

0.1501 

0.1374 

(s.e.) 

(0.0653) 

(0.0667) 

(0.1351) 

(0.0653) 

(0.2570) 

(0.1355) 

0.2 

0.5722 

0.5722 

0.8483 

0.5722 

0.9313 

0.8523 

(s.e.) 

(0.1614) 

(0.1824) 

(0.6283) 

(0.1614) 

(1.4720) 

(0.6290) 

03 

1.1519 

1.1519 

1.7113 

1.1519 

1.8788 

1.7195 

(s.e.) 

(0.2388) 

(0.2945) 

(1.2318) 

(0.2388) 

(2.9542) 

(1.2328) 

A 

9.6187 

— 

8.5587 

9.6187 

3.4632 

3.0614 

(s.e.) 

(1.5569) 

— 

(1.9069) 

(1.5590) 

(2.8867) 

(9.3275) 

P 

— 

9.6187 

0.8117 

2 .1e-12 

0.5548 

211.651 

(s.e.) 

— 

(2.7455) 

(0.2828) 

(0.0455) 

(0.0328) 

(88.5725) 

7 

— 

2.351e-4 

— 

— 

1.2553 

0.8088 

(s.e.) 

— 

1.297e-4 

— 

— 

(0.8749) 

(0.2814) 

6 

— 

— 

— 

— 

0.1462 

— 

(s.e.) 

— 

— 

— 

— 

(1.9357) 

— 

log(£) 

36.670 

36.670 

36.857 

36.670 

36.859 

36.857 

AIC 

-65.340 

-63.340 

-63.714 

-63.340 

-59.717 

-61.714 

AICC 

-64.258 

-61.673 

-62.048 

-61.673 

-56.423 

-59.314 

BIC 

-58.389 

-54.652 

-55.026 

-54.651 

-47.553 

-51.288 

K-S (A x ) 

0.1808 

0.1808 

0.1678 

0.1808 

0.1679 

0.1680 

(p-value) 

(0.1282) 

(0.1282) 

(0.1872) 

(0.1282) 

(0.1869) 

(0.1866) 

K-S (X 2 ) 

0.1410 

0.1411 

0.1289 

0.1410 

0.1290 

0.1291 

(p-value) 

(0.3408) 

(0.3408) 

(0.4499) 

(0.3408) 

(0.4490) 

(0.4484) 

K-S (max(X 1 ,X 2 ) ) 

0.1350 

0.1350 

0.1197 

0.1350 

0.1198 

0.1198 

(p-value) 

(0.3929) 

(0.3929) 

(0.5438) 

(0.3929) 

(0.5428) 

(0.5422) 

LRT 

— 

0.000 

0.374 

0.000 

0.378 

0.374 

(p-value) 

— 

1.0000 

(0.5408) 

1.0000 

0.9447 

0.8294 



6 Conclusions 


In this paper we have introduced the bivariate exponentiated extended Weibull distribution 
whose marginals are exponentiated extended Weibull distributions. We discussed some statis¬ 
tical properties of the new bivariate model. Maximum likelihood estimates of the new class of 
distributions are discussed and we provided the observed Fisher information matrix. Two real 
data sets are used to show the usefulness of the new class. 
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